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DETECTION OF VHE RADIATION FROM THE BL LAC PC 1553+113 WITH THE MAGIC TELESCOPE 
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ABSTRACT 

In 2005 and 2006, the MAGIC telescope has observed very high energy gamma-ray emission from 
the distant BL Lac object PG 1553-1-113. The overall significance of the signal is 8.8 a for 18.8 h 
observation time. The light curve shows no significant fiux variations on a daily time-scale, the 
flux level during 2005 was, however, significantly higher compared to 2006. The differential energy 
spectrum between ~ 90 GeV and 500 GeV is well described by a power law with photon index 
F = 4.2 ± 0.3. The combined 2005 and 2006 energy spectrum provides an upper limit of z = 0.74 on 
the redshift of the object. 

Subject headings: PG 1553-1-113, BL Lac, AGN, VHE gamma-ray, imaging air Cherenkov telescope 



1. INTRODUCTION 

1.1. The BL Lac object PG 1553+113 

The Active Galactic Nucleus (AGN) PG 1553-^113 
was first reported in the Palomar-Green catalogue of UV 
bright objects (Green et al. 1986). It was the only new 
BL Lac object found in the survey and the first BL Lac 
object found in an optical survey. Its spectrum is, typical 
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for BL Lac objects, featureless (Miller & Green 1983) and 
the optical variability strong {mp — 13.2 — 15.0; Miller 
et al. (1988)). The spectral characteristics are close 
to those of X-ray selected BL Lacs (Falomo & Treves 
1990) and it is classified in the literature as intermedi- 
ate BL Lac (Laurent-Muehleisen et al. 1999; Nieppola 
et al. 2006) or high-frequency peaked BL Lac (Giommi 
et al. 1995), as its synchrotron peak frequency lies on the 
borderline of these two groups. 

Despite several attempts, no emission or absorption 
lines have been found in the spectrum of PG 1553-1-113 
(Falomo & Treves 1990). Thus only indirect methods can 
be used to determine the redshift z (e.g. Sbarufatti et al. 
(2005, 2006)). The host galaxy was not resolved in Hub- 
ble Space Telescope (HST) images (Urry et al. 2000), it 
is therefore safe to assume z > 0.25. The observation of 
very high energy (VHE, defined here as E > 100 GeV) 
7-ray emission, on the other hand, may permit to set 
an upper limit on z. The 7-ray absorption in the Extra- 
galactic Background Light (EBL) by means of e~ pair 
production (Stecker et al. 1992; Aharonian et al. 2006a) 
can significantly affect the shape of the observed energy 
spectrum depending on the source redshift. Based on 
present-day EBL models and the observed 7-ray spec- 
trum, one can derive the intrinsic spectrum as a function 
of z. Physical constraints on e.g. the slope of the intrin- 
sic spectrum may then permit to set upper limits on the 
possible redshift (Aharonian et al. 2006b). 

PG 1553+113 belongs to a catalog of X-ray bright ob- 
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jects (Donate ct al. 2005) and, based on its Spectral En- 
ergy Distribution (SED) properties, was one of the most 
promising candidates from a list of VHE 7-ray emitting 
candidates proposed by Costamante & Ghisellini (2002). 
So far, upper limits on the 7-ray emission have been 
reported by the Whipple collaboration (19% Crab flux 
above 390 GcV, de la Calle Perez ct al. (2003)) and Mi- 
lagro (Wilhams 2004). Recently the H.E.S.S. collabora- 
tion has presented evidence for a 7-ray signal at the 4a- 
level (up to 5.3 ct using a low energy threshold analysis) 
above 200 GeV corresponding to about 2% of the Crab 
flux (Aharonian ct al. 2006b) . The energy spcctnmi was 
found to have a steep slope with F = 4.0±0.6stat and an 
upper limit on the redshift of ^; < 0.74 was derived. 

1.2. The MAGIC telescope 

The MAGIC telescope is located on the Canary Is- 
land of La Palma (28.75° N, 17.86° W, at 2225 m 
asL). The telescope comprises a 17 m diameter tessel- 
lated, parabolic mirror with a total surface of 234 m^, a 
light-weight space-frame made from carbon fiber-epoxy 
tubes, and a camera with 576 hemispherical photo- 
multiplier tubes (PMT) with enhanced quantum effi- 
ciency (Pancque et al. 2004). The field of view of the 
camera is 3.5° while the trigger area covers about 2.0° 
in diameter. The fast PMT analog signals are routed 
via optical fibers to the DAQ-system electronics in the 
counting house 80 m away. The signals are digitized by 
dual range 300 MHz FADCs. MAGIC can explore 7-rays 
at energies down to 50 GeV (trigger threshold, depend- 
ing on the zenith angle), critical for the observation of 
medium redshift VHE sources with steeply falling spec- 
tra like PC 1553+113. The MAGIC telescope param- 
eters and performance are described in more detail in 
Baixeras et al. (2004) and Cortina et al. (2005). 
Simultaneous with MAGIC, optical observations were 
performed with the KVA telescope on La Palma, oper- 
ated remotely from Tuorla Observatory. The main in- 
strument is a 60 cm (f/15) Cassegrain telescope equipped 
with a CCD capable of polarimetric measurements. A 
35 cm auxiliary telescope (f/H) is mounted on the same 
RA axis. This telescope is used for BVRI CCD photom- 
etry. 

2. OBSERVATION AND DATA ANALYSIS 

PG 1553+113 was observed with the MAGIC tele- 
scope for 8.9 h in April and May 2005, i.e. at about 
the same time when also H.E.S.S. observed the source, 
and for 19 h from January to April 2006. In addi- 
tion to the observations with MAGIC and the 35 cm 
photometric telescope simultaneous data in X-rays were 
taken with the All-Sky-Monitor on board the RXTE 
satellite. These data are provided on the web at 
http://heasarc.gsfc.nasa.gov/xte_weather/. 
Data taken during non-optimal weather conditions or 
affected by hardware problems were excluded from the 
analysis. Also, only data taken at small zenith angles 
ZA < 30° (corresponding to a low energy threshold 
which is suitable for a steep energy spectrum) were re- 
tained although measurements went up to 53°. After 
these selection cuts, 7.0 h and 11.8 h of good data re- 
mained for 2005 and 2006, respectively. Given the mean 
ZA of ~ 22° 7-ray events above ^ 90 GeV have been 
used for the physics analysis. 



In addition to the so-called on-data from PG 1553+113, 
off-data were taken on a nearby sky position where no 7- 
ray source is expected, but with comparable zenith angle 
distribution and night sky background light conditions. 
The off-data are used to determine the background con- 
tent in the signal region of the on-data. This was done 
by means of a second order polynomial fit (without lin- 
ear term) to the ALPHA distribution of the normalized 
off-data. The normalization was done in the ALPHA 
region between 30° and 90° where no 7-ray events are 
expected. The ALPHA parameter describes the orienta- 
tion of a shower image in the camera with respect to the 
camera center. Air showers which are aligned parallel to 
the telescope axis do have ALPHA values close to zero. In 
total, 14.5 h of off-data (6.5 h from 2005 and 8.0 h from 
2006) have been used for the analysis. Since the two off- 
samples were in good agreement we used the combined 
data to analyze the individual on-data samples. 

The data were analyzed using the standard MAGIC 
analysis programs for calibration, image cleaning, cut op- 
timization and energy reconstruction (Bretz et al. 2005; 
Gaug et al. 2005; Wagner et al. 2005). The primary 
method for discrimination between hadron- and 7-ray- 
induced showers is based on the Random Forest (RF) 
method (Breiman 2001; Bock et al. 2004) which was 
trained on off-data and Monte Carlo (MC) generated 7- 
ray events. The significance of any excess was calculated 
according to Eq. 17 in Li & Ma (1983) where the on to 
off ratio a was derived, taking into account the smaller 
error from the off-data fit. In addition to the cut opti- 
mization, the RF method was also used for the energy 
estimation based on the image parameters of a statisti- 
cally independent MC 7-ray sample. The average energy 
resolution obtained was 24% RMS. All MC data used in 
this analysis were generated using the CORSIKA version 
6.019 (Knapp & Heck 2004; Majumdar et al. 2005). 

3. RESULTS 

Combining the data from 2005 and 2006, a very clear 
signal is seen in the image parameter ALPHA, as shown 
in Fig. 1. Defining the signal region as ALPHA < 12° 

(containing about 90% of the 7-ray events), an excess of 
1032 over 8730 background events yields a total signifi- 
cance of 8.8 (T. The individual results for the years 2005 
and 2006 are listed separately in Table 1. In both years 
the object has been clearly detected with a significance 
> 6ct. 

The 7-ray, X-ray and optical light curve of PG 
1553+113 are shown in Fig. 2. While the optical data 
shows significant short term variability on the 25% level 
the X-ray data is consistent with a constant emission, 
given the weighted mean of 0.15 ± 0.03 counts/s. In 
7-rays there is no evidence for short term variability, 
but a significant change in the flux level from 2005 to 
2006 is found, given a systematic error of the analy- 
sis on the flux level of about 30%. The average in- 
tegral flux between 120 Gc^V and 400 GeV is given as 
F = 10.0±0.23,,tat and F = 3.7 ± 0.08,,taf (with F given 
in miits of lO^'^^ cm^^ s^^) for 2005 and 2006, respec- 
tively. On February 25th, prior to the optical flare, the 
optical polarimetry of the source was measured with the 
KVA 60 cm telescope. The degree of optical linear po- 
larisation was 8.3 ± 0.2% and the polarisation position 
angle was 139.1° ± 0.4°. It should be noted, that since 
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the host galaxy can not be resolved for this object, the 
optical flux should correspond to the emission from the 
AGN core. 

The combined 2005 and 2006 diflferential energy spec- 
trum for PG 1553+113 is shown in Fig. 3. The integral 
fluxes above 200 GeV and the spectral slope coeflicients 
for the different samples are listed in Tab. 1. Effects on 
the spectrum determination introduced by the limited 
energy resolution were corrected by 'unfolding' according 
to Mizobuchi et al. (2005). For comparison, the MAGIC 
Crab energy spectrum and the H.E.S.S. PG 1553+113 
energy spectrum (Aharonian et al. 2006b) are also shown. 
The energy spectrum is well described by a pure power 
law: 



dN 
dE 



(1.8±0.3«tot) 



E 



200 GeV 



-4.2±0.33ta 



(1) 



(in units of 10-i°cm-2 s'^ TeY'^, x /NDF = 1.5/4). 
Compared to the Crab spectrum in the same energy 
range (P = 2.41 ± 0.05; Wagner et al. (2005)) this spec- 
trum is significantly steeper. The spectral slopes of the 
individual years are in good agreement although the flux 
level above 200 GeV is about a factor 3 larger in 2005 
compared to 2006. This is also shown in the light curve. 
The estimated systematic error on the analysis (signal 
extraction, cut efficiencies etc.) is 25% (dark colored 
band in the figure) and 30% on the energy scale (light 
colored band). 

4. DISCUSSION 

The BL Lac object PG 1553+113 has been detected at 
8.8 a with the MAGIC telescope in 18.8 hours of obser- 
vation during 2005 and 2006. This confirms the tentative 
signal seen by H.E.S.S. at a higher energy threshold with 
data taken at about the same time as MAGIC in the 2005 
period (Aharonian et al. 2006b). The source, therefore, 
can now be considered as firmly detected. 
The agreement between the measured H.E.S.S. and 
MAGIC energy spectra of PG 1553+113 in 2005 is rea- 
sonably good. While the spectral slope is consistent 
within errors, the absolute flux above 200 GeV in 2005 
is by a factor 4 larger compared to H.E.S.S. This differ- 
ence may in part be explained by the systematic errors of 
both measurements but also by variations in the flux level 
of the source (the observations with H.E.S.S. were com- 
menced after MAGIC). The observed energy spectrum 
is steeper than that of any other known BL Lac object. 
This may be an indication of a large redshift (z > 0.3), 
but can as well be attributed to intrinsic absorption at 
the AGN or, more naturally, to an inverse Compton peak 
position at lower energies. The spectrum can, however, 
be used to derive an upper limit on the source redshift 
from physical constrains on the intrinsic photon index 
{^int > 1-5) as discussed in Aharonian et al. (2006b). 
Using the lower limit on the evolving EBL density from 
Kneiske et al. (2004) we derived a 2a upper limit on the 
redshift of z < 0.74. The same value was reported by 
Aharonian et al. (2006b) where a slightly different EBL 
model was used. 



The Broad band SED of PG 1553+113 together with 
the results from a model calculation are shown in 
Fig. 4. The VHE data points correspond to the intrin- 
sic spectrum of PG 1553+113 as derived for a redshift 
of z = 0.3. The black points at low energies denote the 
average optical and X-ray flux taken at the same time as 
the MAGIC observations. The gray hatched radio, op- 
tical and X-ray non-simultaneous data were taken from 
Giommi et al. (2002). The solid line shows the result 
of a model fit to the simultaneously recorded data (black 
points) using a homogeneous, one-zone Synchrotron Self- 
Compton (SSC) model as provided by Krawczynski et al. 
(2004). As can be seen from Fig. 4, the 7-ray, X-ray and 
optical data are well described by the model. This is not 
the case for the radio data where intrinsic absorption re- 
quires a much larger emitting volume compared to X-rays 
and 7-rays. Except for a somewhat smaller radius of the 
emitting region, identical model parameters as in Costa- 
mante & Ghisellini (2002) have been used: Doppler fac- 
tor D = 21, magnetic field strength B = 0.7 G, radius of 
the emitting region R = 1.161q'2i " ^^^^ cm, electron en- 
ergy density pe = O-lll^ Jg erg/cm^ slope of the electron 
distribution ae = -2.6 for 8.2 < log {E/eV) < 9.8^"-^^ 

and ae = -3.6 for 9.8^0^0^^ < log{E/eV) < 10.6io;^. 
The limits on some of these parameters indicate the 
change of the SED model parameters when varying the 
assumed redshift from z=0.2 up to z=0.7 (parameters 
without limits were kept constant for all fits). In the 
case of z > 0.56 the SED model can not accurately de- 
scribe the data and, based on the obtained value, a 
redshift of 0.56 is excluded on the 4.5(7 level. For a com- 
parison of the model parameters with those from other 
BL Lacs we refer to Costamante & Ghisellini (2002). 
PG 1553+113 was in a high state in the optical in 
both years showing a strong fiare at the end of March 
2006. The high linear polarization of the optical emis- 
sion (8.3 ± 0.2%) indicates that a sizeable fraction of the 
optical fiux is indeed synchrotron radiation. In 7-rays 
only a significant change in the flux level from 2005 to 
2006 is found while there is no evidence for variability 
in X-rays. As a result, a possible correlation between 
the different energy bands can not be established. A 
possible connection between the 7-ray detection and the 
optical high state can, however, not be excluded. The op- 
tical flare without X-ray or 7-ray counterpart may still 
be explained by external-inverse-Compton (EIC) models 
which predict a time lag of the X-rays and 7-rays with 
respect to the optical emission. 



ACKNOWLEDGMENTS 

We would like to thank the lAC for the excellent work- 
ing conditions on the La Palma Observatory Roque de 
los Muchachos. We are grateful to the ASM/RXTE team 
for their quick- look results. The support of the German 
BMBF and MPG, the Italian INFN and the Spanish CI- 
CYT is gratefully acknowledged. This work was also 
supported by ETH Research Grant TH-34/04-3 and by 
Polish Grant MNil 1P03D01028. 



REFERENCES 



Aharonian, F. et al. 2006a, Nature, 440, 1018 



Aharonian, F. et al. 2006b A&A, 448, L19 



4 



Baixeras, C. et al. 2004, Nucl.Inst.Meth., A 518, 188 

Bock, R. K. et al. 2004, Nucl.Inst.Meth., A 516, 511 

Breiman, L. 2001, Machine Learning, 45, 5 

Bretz, T. et al. 2005, AIP Conference Proceedings, 745, 730 

Costamante L. & Ghisellini, G. 2002, A&A, 384, 56 

Cortina, J. et al. 2005, Proc. 29th ICRC, Pune, 5-359, preprint 

(astro-ph/0508274) 
de la Calle Perez, I. et al. 2003, ApJ, 599, 909 
Donato, D., Sambruna, R. M., & Gliozzi, M. 2005, A&A, 433, 1163 
Falomo, R. & Treves, A. 1990, PASP, 102, 1120 
Gang, M. et al. 2005, Proc. 29th ICRC, Pune, 5-375, preprint 

(astro-ph/0508274) 
Giommi, P., Ansan, S. G., & Micol, A. 1995, A&AS, 109, 267 
Giommi, P., Massaro, E., & Palumbo, G. 2002, Proc. Blazar 

Astrophysics with BeppoSAX and other Observatories, 2002, 63- 

100 

Green, R. F., Schmidt, M., & Licbcrt, J. 1986, ApJS, 61, 305 
Knapp, J. & Heck, D. 2004, EAS Simulation with CORSIKA: A 

Kneiske, T. M.. Bretz, T., Mannheim, K., & Hartmann, D. H. 2004, 

A&A, 413, 807 
Krawczynski, H. et al. 2004, ApJ, 601, 151 

Laurent-Muchlcisen, S. A., Kollgaard, R. I., Feigelson, E. D., 

Brinkmarm, W., & Siebert, J. 1999, ApJ, 525, 127 
Li, T. & Ma, Y. 1983, ApJ 272, 317 



Majumdar, P. et al. 2005, Proc. 29th ICRC, Pune, 5-203, preprint 

(astro-ph/0508274) 
Miller, H. R., Carini, M. T., Gaston, B. J., & Hutter, D. J. 1988, 

Proc. of the lUE Symposium (Grecnbelt), 303 
Miller, H. R. & Green, R. F. 1983, BAAS, 15, 957 
Mizobuchi, S. et al. 2005, Proc. 29th ICRC, Pune, 5-323, preprint 

(astro-ph/0508274) 
Nieppola, E., Tornikoski, M., & Vahaoja, E. 2006, Ai-IA, 445, 441 
Panequc, D., Gcbaucr, H. J., Lorenz, E., & Mirzoyan, R. 2004, 

Nucl.Inst.Meth., A 518 619 
Primack, J., Bullock, J. S., & Sonicrvillc, R. S. 2005, AIP 

Conference Proceedings, 745, 23 
Sbarufatti, B., Treves, A., & Falomo, R. 2005, ApJ, 635, 173 
Sbarufatti, B., Treves, A., Falomo, R., Heidt, J., Kotilainen, J., & 

Scarpa, R. 2006, AJ, 132, 1 
Stecker, F. W., de Jager, O. C, & Salamon, M. H. 1992, ApJ 390, 

49 

Urry, C. M., Scarpa, R., O'Dowd, M., Falomo, R., Pesce, J. E. & 

Treves, A. 2000, ApJ, 532, 816 
Wagner, R. M. et al. 2005, Proc. 29th ICRC, Pune, 4-163, preprint 

(astro-ph/0508244) 
WiUiams, D. 2004, AIP Conference Proceedings, 745, 499 



5 




PG1 553+1 13 

2005/2006 

Observation time: 18.8 h 

Background events: 8730± 66 
Excess Events : 1 032±1 1 9 
Significance: 8.8a 
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Fig. 1. — ALPHA plot for the combined 2005 and 2006 PG 1553+113 data after cuts. The diagram also shows the distribution of the 
(normalized) off-data and a second-order polynomial describing the off-data. 




Fig. 2.— VHE 7-ray (120 GeV - 400 GeV), X-ray (2 keV - 10 keV) and optical light curve (R-Band) of PG 1553+113 in 2005 and 2006. 
The horizontal bars in the top panel correspond to the average flux during 2005 and 2006, respectively. 
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Fig. 3. — Differential energy spectrum of PG 1553+113 as derived from the combined 2005 and 2006 data. Tfic MAGIC Crab energy 
spectrum and the H.E.S.S. PG 1553+113 energy spectrum have been included for comparison. 
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Fig. 4. — Broad band SED of PG 1553+113. The solid lines are the result of a SSC model fit to the black data points using the code 
provided by Krawczynski ct al. (2004) (see text). The gray points comprise non-simultaneous radio, optical and X-ray data from Giommi 
et al. (2002). 
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TABLE 1 

Results from the PG 1553+113 analysis as derived for 2005 and 2006. 



Year 


on time 




No// 


^excess 


on / off 


sigma 


F{E > 200GeV)'' 


photon index 


2005 
2006 
2005+2006 


7.0 h 
11.8 h 
18.8 h 


3944 
5815 
9761 


3501 ± 26 
5228 ± 39 
8730 ± 66 


443 ± 68 
588 ± 86 
1032 ± 119 


0.20 
0.30 
0.49 


6.7 a 
7.0 a 

8.8 a 


2.0±0.6stat±0.6sys 

0.6 ± 0.2stat ± 0.2sys 

1.0 ± OAstat ± O.Ssys 


4.31 ± 0.45 
3.95 ± 0.23 
4.21 + 0.25 



^integral flux in units of 10 cm ^ s ^ 



